Methods for providing generator stator winding ground fault protection

ABSTRACT

The present invention includes a method for providing ground fault protection to a generator. The method includes determining the existing state of the generator and second, determining the updated alarm or trip conditions of a relay connected to the generator in view of the existing state of the generator. This includes estimating a ratio of the magnitude of the 3 rd  harmonic neutral voltage to the magnitude of the 3 rd  harmonic terminal voltage of a generator, along with calculating the residual 3 rd  harmonic voltage from the estimated 3 rd  harmonic ratio. From this, a comparison is made between the energy of the residual voltage with a fractional energy of the 3 rd  harmonic neutral voltage. Based on this comparison, an alarm or trip condition is signaled when the energy of the residual voltage is greater than the fractional energy of the 3 rd  harmonic neutral voltage. The invention further includes adaptive relays and systems utilizing the adaptive relays, along with methods of producing power using the adaptive relays.

CLAIM OF PRIORITY

The present application claims priority to U.S. patent application No. 13/971,433, filed Aug. 20, 2013, now U.S. Pat. No. 9,236,726, issued Jan. 12, 2016, which in turn claims the benefit of provisional application No. 61/691,471 filed on Aug. 21, 2012.

FIELD OF THE INVENTION

The present invention relates to improved methods, devices, and systems for providing ground fault protection for generators and more particularly methods, devices, and systems for providing 100% stator ground fault protection for synchronous generators.

BACKGROUND OF THE INVENTION

Synchronous generators operate within electrical systems to provide uninterrupted power to consumers. A variety of methods are used to protect the generator and thereby the electrical system from internal and external faults. Internal faults in the generator, such as stator winding ground faults, are an important concern to generation utilities worldwide.

One known system and method for 100% stator ground fault protection uses generator 3rd harmonic voltages. In this system and method, operating data from an in-use generator and the larger electrical system are collected and analyzed to set fixed settings into a relay in order to provide stator ground fault protection. These fixed relay settings determine when the relay is to trip or issue an alarm. The fixed relay settings must include extra tolerances (margins) to allow for normal 3rd harmonic voltage variations that occur while operating generators and electrical systems. Moreover, this solution cannot detect stator ground faults in an off-line generator prior to being connected to the electrical system because operating data is required to set the fixed relay settings.

Another known system and method utilizes sub-harmonic voltage injection. However, sub-harmonic voltage injection requires additional equipment to be placed into service, which requires field commissioning. And yet the system and method still require fixed relay settings. Furthermore, the total capacitance-to-ground of the generator stator windings, bus work, and delta-connected transformer windings of the unit transformer windings must be known to ensure the relay settings are correctly determined. Due to higher cost and maintenance of the added equipment, this solution is less attractive compared to the first solution. However, in contrast to the first solution, this solution can detect stator ground faults in an off-line generator prior to being connected to the electrical system.

In addition, known 100% stator ground fault protection systems and methods using 3rd harmonic voltages lack security at low power factors and during unusual power system conditions. Security means the degree of certainty that a protective relay will operate correctly for normal operating conditions of protected equipment. This lack of security causes false tripping or false alarms for the relays and leads to shutdown of the generator, which in turn leads to a costly loss of production and costly testing of the generator. Hence, it is desirable to develop a secure, low-cost stator ground fault protection methods, devices, and systems.

The present invention overcomes one or more of these problems.

SUMMARY OF THE INVENTION

The present invention includes a method for providing ground fault protection to a generator. The method includes determining the existing state of the generator and second, determining the updated alarm or trip conditions of a relay connected to the generator in view of the existing state of the generator. This includes estimating a ratio of the magnitude of the 3^(rd) harmonic neutral voltage to the magnitude of the 3^(rd) harmonic terminal voltage of a generator, along with calculating the residual 3^(rd) harmonic voltage from the estimated 3^(rd) harmonic ratio. From this, a comparison is made between the energy of the residual voltage with a fractional energy of the 3^(rd) harmonic neutral voltage. Based on this comparison, an alarm or trip condition is signaled when the energy of the residual voltage is greater than the fractional energy of the 3^(rd) harmonic neutral voltage. The invention further includes adaptive relays and systems utilizing the adaptive relays, along with methods of producing power using the adaptive relays.

BRIEF DESCRIPTION OF THE FIGURE

FIG. 1 shows a schematic drawing of the system of the present invention.

DETAILED DESCRIPTION

The present invention includes methods of providing 100% stator ground fault protection to generators in general and synchronous generators in particular. The methods allow a relay to adapt to the existing state of the generator to which the relay is connected or to the state of the electrical system to which the generator is connected. The inventive devices (e.g. the adaptive relays) incorporate these methods and the inventive systems incorporated these devices.

Numerous outside factors can influence the state of the generator including the load on the generator, the load on the electrical system, voltage spikes and drops, power swings, etc. Detecting ground faults in the stator winding is difficult in the ever changing state of the generator. So having a better understanding of the instantaneous state of the generator or electrical system is needed. From the state of the generator or electrical system, updated or adaptive alarm or trip conditions can be set on the relay for the generator. In this manner, the generator is protected even while the existing state of the generator and electrical system is taken into account. Recognize that the existing state of the generator takes into account both internal and external influences, such as from the electrical system to which the generator is connected.

The ground fault protection method of the present invention accomplishes these two goals. First, it determines the existing state of the generator and second, it determines the updated alarm or trip conditions of the relay in view of the existing state of the generator. The first goal is accomplished by estimating a ratio of the magnitudes of the 3^(rd) harmonic neutral and terminal voltages and then calculating a residual voltage from the estimated ratio. The second goal is accomplished by comparing the energy of the residual 3^(rd) harmonic voltage with the fractional energy of the 3^(rd) harmonic neutral voltage. Based on this comparison, the existence of an alarm or trip condition is determined and an alarm or trip signal will issue.

Initially, an estimate of the ratio of the magnitudes of the 3^(rd) harmonic voltages of the generator is needed. This estimate is a snap shot of the existing state of the generator. This estimate is in contrast to the direct calculation of the ratio of the magnitudes of the 3^(rd) harmonic voltages. The estimated ratio ({circumflex over (ρ)}) is the ratio of the magnitude of the 3^(rd) harmonic neutral voltage ({tilde over (v)}_(N3)) of the generator to the magnitude of the 3^(rd) harmonic terminal voltage ({tilde over (V)}_(T3)) of the generator, that is,

$\hat{\rho} = {\frac{{\overset{\sim}{V}}_{N\; 3}}{{\overset{\sim}{V}}_{T\; 3}}.}$

The estimated 3^(rd) harmonic ratio is preferred over the directly calculated 3^(rd) harmonic ratio because the estimated 3^(rd) harmonic ratio is less susceptible to measurement and numerical noise and thus leads to a ground fault protection scheme and relay that are more secure.

Any of several iterative techniques are suitable for estimating the 3^(rd) harmonic ratio including weighted least-squares (WLS) and recursive weighted least-squares (RWLS). A RWLS with a ‘forgetting factor’ is preferred (also know as Exponentially-Weighted-Past RWLS or EWP-RWLS). Likewise, use of a Kalman adaptive filter (KAF) may also be used. Each of these techniques is well understood by a skilled artisan and outlined below.

The ratio ({circumflex over (ρ)}) may be estimated using a weighed least squares technique. The ratio may be estimated using the following vector equation:

{circumflex over (ρ)}_(WLS)=[X^(T)WX]⁻¹X^(T)WY (eq. 1), where X is the matrix of {tilde over (V)}_(T3) at time 1 to M, X^(T) is the transposed matrix X, Y is the matrix of {tilde over (V)}_(N3) at time 1 to M, and W is an M by M matrix of the scalar positive weights associated with each error (e(i)). The error (e(i)) is incurred due to noise, disturbance, or simply model mismatches at the i^(th) data point.

The weighted least squares ratio ({circumflex over (ρ)}) may also be estimated using the following scalar equation:

$\begin{matrix} {{{\hat{\rho}}_{WLS} = {\left\lbrack {\sum\limits_{i = 1}^{M}{{w(i)}{{\overset{\sim}{V}}_{T\; 3}^{2}(i)}}} \right\rbrack^{- 1}\left\lbrack {\sum\limits_{i = 1}^{M}{{w(i)}{{\overset{\sim}{V}}_{T\; 3}(i)}{{\overset{\sim}{V}}_{N\; 3}(i)}}} \right\rbrack}},} & \left( {{eq}.\mspace{14mu} 2} \right) \end{matrix}$

where w(i) is scalar positive weights associated with each error (e(i)).

The ratio ({circumflex over (ρ)}) may also be estimated using a recursive weighed least squares technique. The ratio may be estimated using the following algorithm:

$\begin{matrix} {{RWLS}\mspace{14mu} {Algorithm}\text{:}\mspace{14mu} \left\{ {\begin{matrix} {{{{\hat{\rho}}_{WLS}(t)} = {{{\hat{\rho}}_{WLS}\left( {t - 1} \right)} + {{K(t)}{v(t)}}}},} \\ {{{i.c.} = {{\hat{\rho}}_{0} > 0}},} \\ {{{v(t)} = {{{\overset{\sim}{V}}_{N\; 3}(t)} - {{{\overset{\sim}{V}}_{T\; 3}(t)}{{\hat{\rho}}_{WLS}\left( {t - 1} \right)}}}},} \\ {{{K(t)} = {{P(t)}{w(t)}{{\overset{\sim}{V}}_{T\; 3}(t)}}},} \\ {{{P(t)} = {{P\left( {t - 1} \right)} - \frac{{P\left( {t - 1} \right)}{w(t)}{{\overset{\sim}{V}}_{T\; 3}^{2}(t)}{P\left( {t - 1} \right)}}{1 + {{P\left( {t - 1} \right)}{w(t)}{{\overset{\sim}{V}}_{T\; 3}^{2}(t)}}}}},} \\ {{{i.c.} = {\Pi_{0} > 0}},} \\ {{{w(t)} = \frac{1}{{{PF}(t)}^{2} + \delta}},{\delta > 0}} \end{matrix}.} \right.} & \left( {{eq}.\mspace{14mu} 3} \right) \end{matrix}$

Here, PF(t) is the power factor. i.c. is short for initial conditions, {circumflex over (ρ)}₀ is an initial guess for the value of the ratio, Π₀ is the degree of confidence in this guess. A Small Π₀ indicates high confidence while large Π₀ indicates little confidence in the chosen {circumflex over (ρ)}₀. The residual voltage v(t) is the predicted parameter estimation error at time t. It is the difference between {tilde over (V)}_(N3)(t) and its one-step-ahead prediction {tilde over (V)}_(T3)(t) {circumflex over (ρ)}_(WLS)(t−1). Small v(t) implies that {circumflex over (ρ)}_(WLS)(t−1) is good and very little correction is needed. The gain K(t) determines how much v(t) should change the parameter estimate {circumflex over (ρ)}_(WLS)(t−1). The constant δ is chosen to ensure that the weights are large, but finite during low power factor operating points.

The ratio ({circumflex over (ρ)}) may also be estimated using an exponentially weighted past recursive weighed least squares technique. The ratio may be estimated using the following algorithm:

$\begin{matrix} {{EWP}\text{-}{RWLS}\mspace{14mu} {Algoritm}\text{:}\mspace{14mu} \left\{ {\begin{matrix} {{{{\hat{\rho}}_{EWP}(t)} = {{{\hat{\rho}}_{EWP}\left( {t - 1} \right)} + {{K(t)}{v(t)}}}},} \\ {{{i.c.} = {{\hat{\rho}}_{0} > 0}},} \\ {{{v(t)} = {{{\overset{\sim}{V}}_{N\; 3}(t)} - {{{\overset{\sim}{V}}_{T\; 3}(t)}{{\hat{\rho}}_{EWP}\left( {t - 1} \right)}}}},} \\ {{{K(t)} = {{P(t)}{w(t)}{{\overset{\sim}{V}}_{\; {T\; 3}}(t)}}},} \\ {{P(t)} = {\lambda^{- 1}\left\lbrack {{P\left( {t - 1} \right)} - \frac{{P\left( {t - 1} \right)}{w(t)}{{\overset{\sim}{V}}_{T\; 3}^{2}(t)}{P\left( {t - 1} \right)}}{\lambda + {{P\left( {t - 1} \right)}{w(t)}{{\overset{\sim}{V}}_{T\; 3}^{2}(t)}}}} \right\rbrack}} \\ {{{i.c.} = {\Pi_{0} > 0}},} \\ {{{w(t)} = \frac{1}{{{PF}(t)}^{2} + \delta}},{\delta > 0}} \end{matrix}.} \right.} & \left( {{eq}.\mspace{14mu} 4} \right) \end{matrix}$

The variables and constants in this algorithm have the same definitions as in the RWLS algorithm. The difference between the RWLS estimate of the ratio and the EWP-RWLS estimate of the ratio lies λ in the determination of P(t). λ is between 0 and 1 and is the ‘forgetting factor’; that is, λ represents an exponential decaying of earlier and earlier weighted errors (w(t)). When λ=1, then the RWLS and EWP-RWLS equations become equal; that is, there is no decaying of earlier weighted errors.

In addition to the techniques discussed above, the ratio ({circumflex over (ρ)}) may also be estimated using a Kalman adaptive filter (KAF). The adjective adaptive is added to the Kalman filter to indicate the difference between the conventional Kalman filter where the system parameters are time-invariant or perhaps time varying, but known in advance. The ratio may be estimated using the following KAF algorithm:

$\begin{matrix} {{KAF}\mspace{14mu} {Algorithm}\text{:}\mspace{14mu} \left\{ {\begin{matrix} {{{{\hat{\rho}}_{KAF}(t)} = {{{\hat{\rho}}_{KAP}\left( {t - 1} \right)} + {{K(t)}{v(t)}}}},} \\ {{{i.c.} = {{\hat{\rho}}_{0} > 0}},} \\ {{{v(t)} = {{{\overset{\sim}{V}}_{N\; 3}(t)} - {{{\overset{\sim}{V}}_{T\; 3}(t)}{{\hat{\rho}}_{KAF}\left( {t - 1} \right)}}}},} \\ {{{K(t)} = {{P(t)}{w(t)}{{\overset{\sim}{V}}_{T\; 3}(t)}}},} \\ {{P(t)} = {{P\left( {t - 1} \right)} - \frac{{P\left( {t - 1} \right)}{w(t)}{{\overset{\sim}{V}}_{T\; 3}^{2}(t)}{P\left( {t - 1} \right)}}{1 + {{P\left( {t - 1} \right)}{w(t)}{{\overset{\sim}{V}}_{T\; 3}^{2}(t)}}} + N}} \\ {{{i.c.} = {\Pi_{0} > 0}},} \\ {{{w(t)} = \frac{1}{{{PF}(t)}^{2} + \delta}},{\delta > 0}} \end{matrix}.} \right.} & \left( {{eq}.\mspace{11mu} 5} \right) \end{matrix}$

The variables and constants in this algorithm have the same definitions as in the RWLS algorithm. The difference between the KAF estimate of the ratio and the RWLS estimate of the ratio lies in the addition of N in the determination of P(t). The addition of N to the equation prevents P(t) from converging to zero as time t gets large, and is the reason that the adaptation period in the KAF never stops. N plays a similar role to the forgetting factor λ in the EWP-RWLS technique.

With the estimated 3^(rd) harmonic ratio, the residual 3^(rd) harmonic voltage (v) is calculated. The traditional calculation of the residual 3^(rd) harmonic voltage is v(t)={tilde over (V)}_(N3)(t)−{tilde over (V)}_(T3)(t)ρ(t−t_(∝)) (eq. 6), where ρ(t−t₄) is the directly calculated ratio of {tilde over (V)}_(N3) to {tilde over (V)}_(T3). By substituting in the estimated ratio, {tilde over (ρ)}(t−1), the calculation of the residual 3^(rd) harmonic voltage becomes less susceptible to measurement and numerical noises. Thus, the update calculation for the residual 3^(rd) harmonic voltage is: v(t)={tilde over (V)}_(N3)(t)−{tilde over (V)}_(T3)(t){circumflex over (ρ)}(t−1) (eq. 7).

The residual 3^(rd) harmonic voltage represents the estimated balance between the distribution of the 3^(rd) harmonic voltage magnitudes at the neutral versus at the terminal. This estimated balance is an accurate estimate of this distribution throughout the entire operating range of the generator since {circumflex over (ρ)} is constantly updated. And thus, gives an accurate estimate of the existing state of the generator or electrical system.

Now that an understanding of the existing state of the generator (as it may have been influenced by the electrical system) is achieved, the alarm or trip conditions of the relay need to be updated. The alarm or trip conditions are essentially monitoring whether the operate quantity of a relay is larger than the restraint quantity of the relay over a plurality of data points. This is expressed as a comparison of the two quantities.

The traditional comparison is given by |{tilde over (V)}_(N3)(t)−{tilde over (V)}_(T3)(t)ρ(t−t_(∝))|>β{tilde over (V)}_(N3)(t) (eq. 8) and uses the residual 3^(rd) harmonic voltage using the directly calculated ratio of the 3^(rd) harmonic voltages and compares it to a fractional portion of the 3^(rd) harmonic neutral voltage. A better comparison is given by |V_(N3)(t)−{tilde over (V)}_(T3)(t){circumflex over (ρ)}(t−1)|>β{tilde over (V)}_(N3)(t) (eq. 9), where the estimated ratio of 3^(rd) harmonic voltages is used in place of the directly calculated ratio. In both equations, β is a user selected constant between 0 and 1 that gives a fractional portion of the 3^(rd) harmonic neutral voltage. β is selected to modify the sensitivity of comparison calculation. The left side of eq. 9 is the operate quantity and the right side of eq. 9 is the restraint quantity.

When eq. 9 is true, then the alarm or trip condition will exist, and an alarm or trip signal will issue. Stated alternatively, when the operate quantity is greater than the restraint quantity, an alarm or trip signal will be issued by the relay.

This comparison is consistent. That is, for a ground fault at or near the neutral, {tilde over (V)}_(N3) is approximately 0 and {tilde over (V)}_(T3) will be approximately equal to the total 3^(rd) harmonic voltage magnitude. Thus, eq. 9 will collapse down to |{tilde over (V)}_(N3)(t){circumflex over (ρ)}(t−1)|>0 . Thus, the operate quantity will be greater than the restraint quantity, and an alarm or trip will exist. Thus, this permits the detection of ground faults at or near the neutral of the generator. Moreover, this trip or alarm condition will exist as estimated ratio, {circumflex over (ρ)}(t−1) gradually decays to zero.

Also, this comparison permits the detection of ground faults at start up. That is, by setting the starting condition to {circumflex over (ρ)}(0)=1 in eq. 9, the operate quantity will be greater than the restraint quantity.

To improve the security of the ground fault detection scheme, the comparison is carried out multiple times using multiple readings of the 3^(rd) harmonic voltages. These several data points are used to compare energies. That is, energies are calculated from several instantaneous readings of the 3^(rd) harmonic residual and neutral voltages (as shown below). This smoothes out the readings by reducing the impact of any individual instantaneous reading. Thus, the comparison of energies is done with at least three and preferably a minimum of six data points before an alarm or trip condition exists or is signaled.

The operate quantity of the energy is defined as:

$\begin{matrix} {{_{AO}(t)} = \left\{ {\begin{matrix} {0,} & {{{if}\mspace{14mu} 1} \leq t \leq L} \\ {{\sum\limits_{s = {t - L}}^{t}{{v(s)}}^{2}},} & {{{if}\mspace{14mu} t} > L} \end{matrix},,} \right.} & \left( {{eq}.\mspace{11mu} 10} \right) \end{matrix}$

where L is the number of instantaneous readings (i.e. data points) and preferably L is at least six. The restraint quantity is defined as:

$\begin{matrix} {{_{AR}(t)} = \left\{ {\begin{matrix} {{\sum\limits_{s = 1}^{t}{{\overset{\sim}{V}}_{N\; 3}^{2}(s)}},} & {{{if}\mspace{14mu} 1} \leq t \leq L} \\ {{\sum\limits_{s = {t - L}}^{t}{{\overset{\sim}{V}}_{N\; 3}^{2}(s)}},} & {{{if}\mspace{14mu} t} > L} \end{matrix}.} \right.} & \left( {{eq}.\mspace{14mu} 11} \right) \end{matrix}$

These two quantities, representing energies, are compared:

AO(t)>β

AR(t) (eq. 12), and an alarm or trip is signaled if this equation is satisfied for any t>L. Note that the window of the L samples is needed to ensure that the adaptive ratio estimation algorithm has properly learned the current operating conditions of the generator or electrical system. Moreover, depending on the timing of a fault which occurs after the learning period is completed, an alarm or trip signal is issued, if over a finite window of L samples, the energy of the residual signal v(t) is strictly greater than the fraction of the energy of {tilde over (V)}_(N3)(t).

The present invention also includes devices, called adaptive relays that carry out the invention method. The adaptive relay includes a processor, which in turn is programmable. Programming of the processor permits the alarm or trip conditions of the relay to be modified on a real time basis; that is, the trip or alarm conditions for the relay can be continuously adjusted to meet the existing state of the generator and the electrical system. The continuously adjusted relay settings provide a secure and dependable protection to the generator. Exemplary programmable relays include SEL-300G and SEL-700G relays.

Essentially, the processor in the relay is programmed to carry out the methods discussed above; namely, determining the existing state of the generator and electric system and then updating the alarm or trip conditions of the relay.

In addition to the adaptive relay, the present invention also relates to systems that incorporate the adaptive relay, as shown in FIG. 1; in particular, those systems 10 in which one or more adaptive relays 12 are electrically connected to a generator 14.

The method of operating the generator includes producing power with a generator having an adaptive relay electrically connected to the generator.

It will be further appreciated that functions or structures of a plurality of components or steps may be combined into a single component or step, or the functions or structures of one-step or component may be split among plural steps or components. The present invention contemplates all of these combinations. Unless stated otherwise, dimensions and geometries of the various structures depicted herein are not intended to be restrictive of the invention, and other dimensions or geometries are possible. Plural structural components or steps can be provided by a single integrated structure or step. Alternatively, a single integrated structure or step might be divided into separate plural components or steps. In addition, while a feature of the present invention may have been described in the context of only one of the illustrated embodiments, such feature may be combined with one or more other features of other embodiments, for any given application. It will also be appreciated from the above that the fabrication of the unique structures herein and the operation thereof also constitute methods in accordance with the present invention. The present invention also encompasses intermediate and end products resulting from the practice of the methods herein. The use of “or” should be understood to include by the conjunctive and disjunctive definitions of the word. The use of “comprising” or “including” also contemplates embodiments that “consist essentially of” or “consist of” the recited feature.

The explanations and illustrations presented herein are intended to acquaint others skilled in the art with the invention, its principles, and its practical application. Those skilled in the art may adapt and apply the invention in its numerous forms, as may be best suited to the requirements of a particular use. Accordingly, the specific embodiments of the present invention as set forth are not intended as being exhaustive or limiting of the invention. The scope of the invention should, therefore, be determined not with reference to the above description, but should instead be determined with reference to the appended claims, along with the full scope of equivalents to which such claims are entitled. The disclosures of all articles and references, including patent applications and publications, are incorporated by reference for all purposes. 

What is claimed is:
 1. A relay, comprising: a processor capable of being programmed with alarm or trip conditions for a relay; and a program for the processor that establishes the conditions for the relay based on an existing state of a generator or electrical system.
 2. The relay of claim 1 wherein the program determines the existing state of the generator or electrical system by estimating a ratio of the magnitude of the 3^(rd) harmonic neutral voltage to the magnitude of the 3^(rd) harmonic terminal voltage.
 3. The relay of claim 1, wherein the program determines the existing state of the generator or electrical system by estimating a ratio of the phasor of the 3^(rd) harmonic neutral voltage to the phasor of the 3^(rd) harmonic terminal voltage.
 4. The relay of claim 2 wherein the program calculates a residual 3^(rd) harmonic voltage.
 5. The relay of claim 4 wherein the program establishes modified characteristic of a generator or electrical system relay based on the calculated energy of the residual 3^(rd) harmonic voltage over at least six data points as compared to a fraction of the calculated energy of the 3^(rd) harmonic neutral voltage over the same data points.
 6. The relay of claim 5 wherein the program establishes the conditions for the relay in real time.
 7. The relay of claim 2 wherein estimating the ratio comprises an iterative technique to estimate the ratio of the magnitude of the 3rd harmonic neutral voltage and the magnitude of the 3rd harmonic terminal voltage.
 8. The relay of claim 7 wherein the iterative technique is selected from recursive weighted least-squares (RWLS), exponentially-weighted-past RWLS (EWP-RWLS) or a Kalman adaptive filter (KAF).
 9. An article of manufacture, comprising: a programmable relay, wherein the relay signals an alarm or trip condition when an energy of a residual 3rd harmonic voltage n the existing state of a generator or electrical system including a generator is greater than a fractional energy of a 3rd harmonic neutral voltage in the existing state of a generator or electrical system including a generator.
 10. The article of claim 7 wherein the alarm or trip condition of the relay is set in real time.
 11. The article of claim 7 wherein the residual 3rd harmonic voltage is calculated from an estimated 3rd harmonic ratio.
 12. The article of claim 9 wherein the estimated 3rd harmonic ratio is estimated by comparing a magnitude of the 3rd harmonic neutral voltage to a magnitude of the 3rd harmonic terminal voltage of the generator.
 13. The article of claim 12 wherein estimating the ratio comprises an iterative technique to estimate the ratio of the magnitude of the 3rd harmonic neutral voltage and the magnitude of the 3rd harmonic terminal voltage.
 14. The article of claim 13 wherein the iterative technique is selected from recursive weighted least-squares (RWLS), exponentially-weighted-past RWLS (EWP-RWLS) or a Kalman adaptive filter (KAF).
 15. The article of claim 7 wherein the estimated 3rd harmonic ratio is calculated by using an adaptive algorithm wherein the generator power factor, the magnitude or phasor of the 3rd harmonic neutral voltage, and the magnitude or phasor of the 3rd harmonic terminal voltage are used to estimate the 3rd harmonic ratio. 